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(N ■ ABSTRACT 
(N , 

, We study in detail the kinematics of semileptonic M13 (M —* M'lv) decay, 

considering only the new models preserving all the gauge symmetries of the 
Standard Model (SM) , and assuming that neutrinos are massless and left-handed 
as in the SM. And we present a brief review of the recent theoretical progress on 
model independent constraints on the form factors of M13 decays using dispersion 
relations. Finally we review a new parton model approach for semileptonic 
Qh! B — > D(D*)lv decays, by extending the inclusive parton model and by combining 

with the results of the heavy quark effective theory. We also obtain the slope of 
the Isgur-Wise function. 



^ ■ 1. Introduction 



In exclusive weak decay processes of hadrons, the effects of strong interaction 
are encoded in hadronic form factors. These decay form factors are Lorentz invari- 
ant functions which depend on the momentum transfer q 2 , and their behaviors with 
varying q 2 are dominated by non-perturbative effects of QCD. 

Over the past few years, a great progress has been achieved in our understanding 
of the exclusive semileptonic decays of heavy flavors to heavy flavors^. In the limit 
where the mass of the heavy quark is taken to infinity, its strong interactions become 
independent of its mass and spin, and depend only on its velocity. This provides a 
new SU(2Nf) spin-flavor symmetry, which is not manifest in the theory of QCD. 
However, this new symmetry has been made explicit in a framework of the heavy 
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quark effective theory (HQET)i. In practice, the HQET and this new symmetry 
relate all the hadronic matrix elements of B — ► D and B — ► D* semileptonic decays, 
and all the form factors can be reduced to a single universal function, the so-called 
Isgur-Wise (IW) function!'!, which represents the common non-perturbative dynamics 
of weak decays of heavy mesons. However, the HQET cannot predict the values of the 
IW function over the whole q 2 range, though the normalization of the IW function is 
precisely known in the zero recoil limit. Hence the extrapolation of q 2 dependences 
of the IW function and of all form factors is still model dependent and the source 
of uncertainties in any theoretical model. Therefore, it is strongly recommended to 
determine hadronic form factors of Mi% (M — > M'lv) decay more reliably, when we 
think of their importance in theoretical and experimental analyses. 

Although the Standard Model (SM) has been an enormous success in explaining 
all the known experimental data, it is generally anticipated that there is new physics 
in higher energy regions. One of the long-shot efforts to exploring such new physics 
could be searches for CP violation (or T violation) outside of the Cabibbo-Kobayashi- 
Maskawa (CKM) paradigm!. It is well known! that measuring a component of muon 
polarization normal to the decay plane in K^ 3 (K + — > ■n Q ji + u) decay would signal T 
violation. It is also expected that the CKM phase does not induce the perpendicular 
muon polarization in decay. Therefore, measurements! of these polarizations 
could be clear signatures of physics beyond the SM. 

In Section 2, we study in detail the kinematics of semileptonic M; 3 decay, con- 
sidering only the new models preserving all the gauge symmetries of the SM, and 
assuming that neutrinos are massless and left-handed as in the SM. And in Section 
3, we present a brief review of the recent theoretical progress on model independent 
constraints on the form factors of M/3 decays using dispersion relations. Finally in 
Section 4, we review a new parton model approach for semileptonic B — > D(D*)lu 
decays, by extending the inclusive parton model and by combining with the results 
of the HQET. We also obtain the slope of IW function. 



2. Kinematics of semileptonic M Vi decay 

The semileptonic hadron decays are very important in (i) determining the CKM 
matrix within the SM, and (ii) investigating new effects from physics beyond the SM. 
In this report, considering only the new models preserving all the gauge symmetries 
of the SM, and assuming that neutrinos are massless and left-handed as in the SM, we 
concentrate on the semileptonic decay of a pseudo-scalar meson M to a pseudo-scalar 
meson M' in the M rest frame 

M(p) ^ M'(p')+l(k) + u(k'), (1) 

where all the four-momenta are given in the parentheses. In the SM the semileptonic 
decay process (|l|) proceeds at quark level through a W boson exchange, while the 
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decay process may proceed through charged scalar/tensor exchanges in new physics 
beyond the SM as well. 

The decay amplitude for the process M — > M'lv can be parametrized in the 
factorized form as 



M oc GpV 



(2) 



where Gp is the Fermi decay constant and V is a CKM matrix element. While the 
leptonic current, Z/, is completely determined, the determination of the hadronic cur- 
rent, H^, is limited due to the complicated strong interactions. However, translation 
and Lorentz invariances force if M to be of the form for a scalar or vector exchange: 



(p + p'rf+(t) + (p-pTf-(t), 
P ' q q")f(t) + ^d(t), 



V 



Q- 



(3) 



where q = p — p' = k + k' and t = q 2 . In Eq. (|3|) we have introduced two different 
parametrizations, which are related as follows 



f{t) = 2/+(f), d(t) = (m 2 - m' 2 )f + {t) + t/_(f), 



(4) 



where m(m') is the mass of the pseudo-scalar meson M(M'). Clearly, the so-called 
scalar form factor d(t) vanishes in the limit where the hadronic current is conserved. 
For our convenience, the latter form-factor set (f(t),d(t)) will be employed in the 
followings. 

It is now straightforward to calculate the absolute square of the decay amplitude 
(fj), which can be cast into the form 



|X| 2 ocG||^ 12 | 2 



T f \f\ 2 + T d \d\ 2 + T fd K(fd*) + A fd l(fd*) 



(5) 



In general, the distributions, Tj (i = f,d,fd) and Af d are functions of the four- 
momenta and the lepton polarization vector, which can be explicitly determined. 
First of all, we note that, when the summation over polarization of the lepton is 
taken, the distribution, Af d , vanishes and the other distributions are given by 



Tf 

T d 
Tfd 



8t 2 



— \t(t - m]) Ut + m 2 - m' 2 ) 2 - Am 2 t 



(t + mf){t + m 2 - m' 2 ) - Atx 



1 

4t2 



(t + m 2 )(t + m 2 - m' 2 ) - Atx 



(6) 
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where x = mEi with the lepton energy Ei in the M rest frame. Secondly, we can 
see that three terms \f\ 2 , \d\ 2 , and 1Z(fd*), are independently separable due to the 
different dependence of the three distributions on the variable x. Of course, in this 
case, the lepton mass mi should not be so small in order for the |<i| 2 term to be 
measurable. On the other hand, the differential decay width is given in terms of t 
and x by 

1 \M\ 2 

where m 2 < t < (m — m') 2 and the scatter plot in t and x, which is called a Dalitz 
plot, has uniform phase space density. 

With the lepton polarization into account, we find that the distribution Afj in 
the M rest frame is nothing but the so-called triple vector correlation 



2m;m 



s r (kx p) 



(8) 



where Si is the lepton polarization vector. 

T-invariance requires that the relative phase of / and d should be the same or, in 
other words, I(fd*) should be zero. Conversely, any non-zero value of I(fd*) would 
signal violation of T-violationi. It should be noted that the T-violation effect can not 
be measured without the triple vector correlation (H). Moreover, the T-violation effect 
is proportional to the mass of the lepton and requires the measurement of the the 
lepton polarization. While the polarization of the electron is almost not measurable, 
the polarization of the muon or tau leptons can be easily determined via the angular 
or energy distributions of their decay products. In light of these features, it is clear 
that in the K meson case is by far more sensitive than and in the B meson 
case the r lepton mode is most sensitive to T(fd*). 



3. Model-independent constraints on the form factors 

The form factors at specific kinematic points can be approximately determined by 
the heavy quark symmetries for heavy meson decays and by the chiral symmetries for 
light meson decays. However, the event rate vanishes at the kinematic point so that 
the extrapolation of the form factor to other kinematical values QCDis required to 
determine the CKM matrix with good precision. Moreover, the form factor extrapola- 
tion is unavoidably accompanied by an uncertainty due to the choice of parametriza- 
tion. Estimates of this uncertainty obtained by varying parametrizations suffer the 
same ambiguity as well. Recently, the possibility of obtaining model-independent 
bounds on the form factors f(t) and d(t) using dispersion relations has revived much 
interestBJl. In this section we present a brief review of the model-independent con- 
straints on meson form factors using dispersion relations. 
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The pointi'00 is that the form factors fit) and d[t) of the semileptonic decay 
process M —>■ M'lvi can be constrained by the knowledge of the two-point function 

i f d 4 xe iq - x (0\T(V^x)V uj (0)\0) = -{g^t - <f<f)II T (t) + g^U L (t), (9) 



in the deep Euclidean region, where V^(x) denotes the vector current responsible for 
the transition M — » M'. In order to render both sides of the relation (^) finite we 
employ the once-subtracted dispersion relations 

X**® ) = -W~\*~* = « Jo (10) 

where the spectral functions ImIIr,x(^) are defined by the relation 
-{g^t - q»q v )\mR T (t) + ^ImII L (0 

= ^/Edpr(27r) 4 %-p r )(0|V(0)|r)(r|^t ( o)|o), (11) 

with the summation over all possible hadron states V of equal flavour quantum num- 
bers and with an integral over the phase space of allowed intermediate states. The 
hadron states V might be resonances or open MM' states in a specific isospin chan- 
nel determined by the quantum numbers of the relevant vector current V\ In some 
cases the open MM' state is precisely the lowest hadronic state contributing to the 
absorptive amplitude. Up to an overall Clebsh-Gordon coefficient T], the function 
(0 1 (0)| MM') is the same analytic function as {M'\V^{Q)\M) which appears in the 
decay M -> M'lv h i.e., 

(0|V(0)|MM'> = J - f(t) + ^d(t) J, (12) 

where the regime of the variable t is (m + m') 2 <t<oo. 

A judicious choice of the indices \i and v makes the spectral functions a sum of 
positive definite terms, so one can obtain strict inequalities by concentrating on the 
term with intermediate states of MM' pairs. In particular, the intermediate open 
state MM' gives the contribution 



1 '\d(t)\ 2 + ( 1 _t_\( 1 _ m mi* 



2 

imn T (t) > JL-yJ(t-t_)(t-t+ 



t 2 V t J V t 



,2 



ImH L (t) > -?l-J(t-t^(t-t + )\d(t)\ 2 , (13) 

for t > t + with t + = (m + m') 2 and t_ = (m — m') 2 . Clearly, the inequalities become 
tighter with more possible contributions to the right side, which are always positive. 
Eqs. ( [TDD an d (0) a U ow us to obtain the inequalities on Xt,l(Q 2 ), e.g., 

V 2 r°° \/(* -*-)(*-*+) 
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The spectral function Xt,l{Q 2 ) can be computed reliably from perturbative for Q 2 far 
from the resonance region and their lowest-order expressions from a quark one-loop 
diagram are 

i , 2x 2 (l - x) 2 



x(l — x)Q 2 + m 2 x + m 2 ,{l 



x 



q~ ' q' 

j>. 3 f 1 x{l — x)[m q x — m q i{l — x)] 



Xl(Q ) = 7^(m, - rn q >) / dx— — — — ■ -, (15) 

47r Jo x(l — x)Q 2 + m l q x + m q ,(l — x) 

where the quark q (q') is inside M (M'). 




Fig. 1. Mapping of the cut i-plane into the unit disc by the conformal transformation. 

Using knowledge of the analytic structure of the form factors plus the bounds ( |l3|) 
one can derive model-independent bounds on the form factors f(t) and d(t) in the 
physical region of semileptonic decay, mf <t< (m — m') 2 . To this end let us map the 
complex t plane onto the unit disk \z\ < 1 by J {t + — £)/ (t+ + Q 2 ) — (1 + z)/(l — z) 
(For the mapping, see Fig. 1), and then we can express J{Q 2 ) as the norm squared 
of an analytic function, 

J(Q 2 ) = = J_ / 27r d%(e*V(e* e )| 2 , Hz) = (p{z)d(z), (16) 

Z7T JO 

on the unit disc of the complex z-plane with an explicitly determined function <p(z). 
With the inner product of functions defined on the unit disc as 

1 r 27r 

M*)Mz)) = — d9gt(e i6 )g 2 (e id ), (17) 

Z7T JO 
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we find that, for gi{z) = 1/(1 — X*z) and #2(2) = h(z), (gi(z)\h(z)) = d(X)ip(X), 
d can be determined" with an appropriate A at any point of the physical region 
mf <t<t- for the semileptonic decay M — *> M'lv\. Furthermore, the 2x2 positive 
semidefinite matrix of Aij = (gi(z)\gj(z)), satisfying det(Aij) > 0, provides us with 
the model- independent bound on d(X) 



|rf(A)l ^ - OTRA)F " (1-|A| 2 )|^(A)P- (18) 



The bounds on the form factor f(t) can be obtained through the same approach as 
well. 

To summarize, the main aspects of the model-independent perturbative QCD 
method are as follows. 

• The monemtum transfer squared Q 2 should be far from the resonance region for 
reliable perturbative QCD calculations. For example, for the heavy B meson 
Q 2 = can be takeni because of large resonance mass compared to the QCD 
scale. However, for the light K or n mesons, the value of Q 2 should be large!. 

• The spectral functions can be further saturated with more contributions in- 
cluded. However, the inclusion of resonancesBl can cause uncertainties due to 
no precise information on their decay constants and masses. Of course, res- 
onances with much larger than t + can be safely ignored. In this light, the 
K — > % mo del is a very good system because of no resonances for t < t + . The 
B — ► 7r mode involvesQ one resonance B*, while the B — » D mode involves a 
lot of resonances, which should be included for reliable predictions on the form 
factors. 

• It is naturally expected that any theoretical or experimental inputsB! of the 
form factors f(t) and d(t) within the physical region should give much stronger 
constraints on the form factors. 

As briefly reviewed so far, the pertuarbatively calculable two-point functions in 
the deep Euclidean region provide us with model-independent bounds on the form 
factors of semileptonic meson decays and the constraints can get much more stringent 
with more theoretical and experimental determinations of the form factors at specific 
kinematic points. A systematic investigation along this line is in progress^. 

4. Parton model approach for semileptonic B — > D(D*)lu 

In this part we developed the parton model approach for exclusive semileptonic 
B decays to D, D*, and predicted the q 2 dependences of all form factors. For the 
details of this exclusive parton model approach, see Ref. Kim et al.El. Previously 
the parton model approach has been established to describe inclusive semileptonic B 
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decaysEl, and found to give excellent agreements with experiments for electron energy 
spectrum at all energies. 

While many attemptsEl'0 describing exclusive B decays often take the pole- 
dominance ansatze as behaviors of form factors with varying q 2 , in our approach they 
are derived by the kinematical relations between initial b quark and final c quark. 
According to the Wirbel et al. mo deli, which is one of the most popular models to 
describe exclusive decays of B mesons, the hadronic form factors are related to the 
meson wavef unctions' overlap-integral in the infinite momentum frame, but in our 
model they are determined by integral of the fragmentation functions, which are at 
least experimentally measurable. 

We now develop the parton model approach for exclusive semileptonic decays of B 
meson by extending the inclusive parton model, and by combining with the results of 
the HQET. Theoretical formulation of this approach is, in a sense, closely related to 
Drell-Yan process, while the parton model of inclusive B decays is motivated by deep 
inelastic scattering process. And the bound state effects of exclusive B decays are 
encoded into the hadronic distribution functions of partons inside an initial B meson 
and of partons of a final state resonance hadron. Then, the Lorentz invariant hadronic 
decay width can be obtained using the hadronic distribution functions /d,_b(x), 

E b ■ dT(B -> D(D*)ev) = J dx J dy f B (x) E b ■ dT(b -> ceu) f D (y) . (19) 

The first integral represents the effects of motion of b quark within B meson and the 
second integral those of c quark within D meson. The variables x and y are fractions 
of momenta of partons to momenta of mesons, 

Pb = xp B , Pc = VPD , (20) 

in the infinite momentum frame. The functions /b(x) and fo{y) are the distribution 
function of b quark inside B meson, and the fragmentation function of c quark to D 
meson respectively. Since the momentum fractions and the distribution functions are 
all defined in the infinite momentum frame, we have to consider the colli derLorentz 
invariant quantity, E ■ dT, to use at any other frame. 

The distribution function can be identified with the fragmentation function for a 
fast moving 6-quark to hadronize into a B meson in the infinite momentum frame. 
In general the distribution and fragmentation functions of a heavy quark (Q = t, b, c) 
in a heavy meson (Qq), which are closely related by a time reversal transformation, 
are of very similar functional forms, and peak both at large value of x. Brodsky et 
al.§ have calculated the distribution function of a heavy quark, which has the same 
form as Peterson's fragmentation function0. Therefore, here we follow the previous 
wor ksS to use the Peterson's fragmentation function for both distributions, Jb{x) 
and fo{y)- It has the functional form: 

f Q (z) = N Q z- 1 (l-I-^) 2 , (21) 
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where Nq is a normalization constant, and Q denotes b or c quark. This functional 
form is not purely ad hoc, but motivated by general theoretical arguments, that 
the transition amplitude for a fast moving heavy quark Q to fragment into a heavy 
meson (Qq) is proportional to the inverse of the energy transfer AE^ 1 . When we 
use this function, we have the advantage that the parameter eg of this function has 
been determined from high energy experiments. Thus we can replace uncertain wave 
functions of heavy mesons by the experimentally measurable fragmentation functions. 

From Lorentz invariance we write the matrix element of the decay B — > Dev in 
the form 

<D\J fl \B>=f + (q 2 )(p B +p D ) IM + f4q 2 )(p B -p D ) IM , (22) 

and in terms of the HQET 

< D{v r )\J,\B{v) >= • v')(v + v% + £-(v v')(v - v%) . (23) 

Due to the conservation of leptonic currents, the form factors multiplicated by q^ do 
not contribute. Then the hadronic tensor is given by 

H„ v = < D\J^\B >< D\J U \B >* 

= 2\f + {q 2 )\ 2 {p B ^p Du +p Bu p D ^) , (24) 

and can be expressed by the IW function, 

H^ = R- l Uv-v')\ 2 {p B ^p Du +p Bv p D ^)(^ + 0{^-)^ , (25) 

where 

_ 2^/m B m D 
m B + m D 

By comparing the matrix element of the decay B — > Dev in the parton model 
approach and the matrix element in terms of the HQET, we can define the function 
F(q 2 ) as 

F(q 2 ) = Jdx f B (x)f D (y(x,q 2 )) xy 3 (x,q 2 ) . (26) 

For given q 2 in our parton picture, the function F(q 2 ) measures the weighted transi- 
tion amplitude, which is explicitly given by the overlap integral of distribution func- 
tions of initial and final state hadrons. Comparing the hadronic tensor of HQET and 
that of parton model approach, the IW function is calculated 

Uv-v')\ 2 (l + 0(^)\ =AN -R-F(v-v') . (27) 
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Finally the q 2 spectrum is given by 

- G iM^ (« - -* + ^ - • <* 

We can also obtain the decay spectrum for B — > D*lv 

dT(B^D*eu) G 2 F \V cb \ 2 , 2 , f , 2 2 , 2 , 2 . 2 2 xi/2 



2 tjr / 2\ // 2 2 , 2\2 2 



+ \rn 2 B W 2 (q 2 ) (m% - m 2 D , + q 2 ) + 3m 2 B W 3 {q 2 



,(29) 



where 

W^q 2 ) = -iVxfl 



2g 2 



^ 2 (g 2 ) = iV^ml-m^ + g 2 ) ( 1 - ^ ^ ] - 2iV 3 g 2 ^ 1 ' /_ 



{m B + m D ,) 2 ) " \ {m B + m D ,) 2 



W 3 (q 2 ) = -N im 2 B q 2 {l 2<r 



(m B + m D * 



\2 



+ N 3 q 2 {m 2 B -m 2 D ,+q 2 )(\ 



+ N 2 q 2 (m 2 B + m 2 D «-q 2 )(l- 7 ^ -) , (30) 



(m B + m D *) 2 / 
2q 2 

v (m B + m D *) 2 / 

and J-'iq 2 ) is defined in (PUp. 

The result is plotted in Fig. 2, also compared with the recent CLEO data0. The 
thick solid line is our model prediction with the parameters (e& =0.004, e c =0.04) 
for the heavy quark fragmentation functions^, the thin solid line the Wirbel et al. 
model prediction0, and the dotted line the Korner et al. model predictionEl. We also 
obtain the values of the slope parameter p within the parton model framework, 

p 2 = 0.582 - 0.896 , 

which are compatible with the average value measured by experiments^, 

p 2 = 0.87 ±0.12 . 

For more details of predictions from this exclusive parton model approach, see Ref. 
Kim et alM 



10 



0.2 

0.18 - 

0.16 - 

0.14 - 

0.12 - 



l dr 

r dg 2 



t 1 1 1 r 



i i 



<> 



<> 



<> 




Fig. 2. q 2 spectrum in B — > D*ev decays. 
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